X(3872) as a hybrid state of the charmonium and the hadronic molecule 
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In order to understand the structure of the X(3872), the cc charmonium core state coupling to 
the D°D*° and D + D*~ molecular states are studied. The strengths of the couplings between the 
charmonium state and the hadronic molecular states are determined so as to reproduce the observed 
mass of the the X(3872). The DD* attraction is taken to be consistent with the observed Z(,(10610) 
and Zt(10650) masses. The isospin symmetry breaking is introduced by the mass differences of the 
neutral and the charged D mesons. The obtained structure of the X(3872) is that about 6% is cc 
charmonium, 69% is the isoscalar DD* molecule and 26% is the isovector DD* molecule, which 
explains many of the observed properties of the X(3872), such as the isospin symmetry breaking, 
the production rate in the pp collision, a lack of the existence of the Xci (2-P) peak predicted by the 
quark model, and the absence of the charged X. The structure of the X(3872) we have obtained is 
not just a D° D*° hadronic molecule but the charmonium-hadronic molecule hybrid state. 

PACS numbers: 14.40.Rt, 12.40.Yx, 13.25. Gv 
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I. INTRODUCTION 

The X(3872) state was first observed in 2003 by Belle 
in — > J/i/jtt + n~ K ± Jl| and was confirmed by CDF 
0, DO (3 and BABAR [J collaborations. The observed 
masses of the V(3872) in the J / 'ip ir + n~ channel from 
the recent measurements of the charged and the neutral 
B decays are (3871.4 ± 0.6 ± 0.1) MeV and (3868.7 ± 
1.5 ± 0.4) MeV, respectively Those from the pp and 
the pp collisions are (3871.61 ± 0.16 ± 0.19) MeV H and 
(3871.95 ± 0.48 ± 0.12) MeV 0). The average mass given 
by the particle data group in 2012 Q is (3871.68 ± 0.17) 
MeV, which is 0.16 MeV below the D°D*° threshold. 
The full width is less than 1.2 MeV. 

As for the spin-parity quantum numbers of the 
V(3872), the angular distributions and correlations of 
the Tr + ir~J/ip final state have been studied by CDF [9( 
and they concluded that the pion pairs originate from 
p° mesons and that the favored quantum numbers of the 
A"(3872) are J PC = 1 ++ and 2 '". Recent analyses sup- 
port J PC = 1++ interpretation [loj]. BABAR has found 
the evidence of the radiative decays of A" (3872) — > jJ/ip 
with 3.4-3.6 a significance [H|,[ll], which implies that the 
C-parity of A" (3872) is positive. Though we assume that 
A(3872) is J PC = 1 ++ , whether the quantum number is 
1 ++ or 2 h is still an issue of the discussion and more 
experimental data are certainly necessary. 

Since the first observation of the A" (3872), it has re- 
ceived much attention because its features are difficult 
to explain if a simple cc bound state of the quark po- 
tential model is assumed [HI]. A" (3872) is one of the 



promising candidates of the exotic states reviewed in 
Ref. [Ull - tlq . Many kinds of structures have been sug- 
gested for the A (3872) from the theoretical side, such as 
a tetraquark structure [T^-HOl , a D°D*° molecule HH- 
I2H ] and a charmonium-molecule hybrid [29l - [3l1 |. We also 
employ this hybrid picture and argue that that is most 
appropriate. 

One of the important properties of the A (3872) is its 
isospin structure. The branching fractions measured by 
Belle 32] is 



Br{X tt+tt-tt J/ip) 
Br(X tt+tt-J/i/j) 



1.0 ±0.4 ±0.3, (1) 



and (0.8 ±0.3) by BABAR [Hj]. Here the two-pion mode 
originates from the isovector p meson while the three- 
pion mode comes from the isoscalar u> meson. So, the eq. 
([1]) indicates strong isospin violation. M. Suzuki has es- 
timated the kinematical suppression factor including the 
difference of the vector meson d ecay width and obtained 
the production amplitude ratio [34| using Belle's value 



A(pJ/i>) 



A{ujJ/yb) 



= 0.27 ±0.02. 



(2) 
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Usual size of the isospin symmetry breaking is at most a 
few %. It is interesting to know what is the origin of this 
strong isospin symmetry breaking. In 35J, this problem 
has been studied by using the chiral unitary model and 
the effect of the p-uj mixing has been discussed in [36| . It 
was reported that both of the approaches can explain the 
observed ratio given in Eq. (JlJ, but at present, consensus 
on the mechanism of the large isospin symmetry breaking 
has yet to be reached. We will show in this work that the 
mass difference of the D°D*° and the D+D*- thresholds 
gives enough amount of the isospin violation to explain 
the experiments. 
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The hadronic decays of the A (3872) are investigated 
in [37i - f45l | and the radiative decays in p6l - l53 |. which de- 
pend on rather detailed of the wave function. They are 
beyond the scope of this paper and we will discuss the 
matter elsewhere. Let us mention briefly, however, that 
the production processes have been studied in [55l - l58l |. 
and it seems that a pure molecule picture cannot explain 
the production rate of the A (3872) in the pp collision well 
[57| . There, the production rate of the A (3872) is about 
1/20 of the rate of -0(25"), which suggests that A(3872) 
has to have, very roughly speaking, the order of 5% of 
the cc component. 

Recently, Z b (10610) ± '° and Z b (10650) ± '° (J p =l+) 
resonances have been found in the T(5S) decay 
to T(riSV+7r- (n=l,2,3) and h b (mP)n + n~ (m=l,2) 
reactions [59]. The masses of these resonances are just 
above to the BB* and the B*B* thresholds, respectively; 
the main components are considered to be the B^B* 
two-meson states. This means that there exist almost- 
zero-energy bound states (or resonances) in both of the 
DD* and the B^B* systems. In order to make such 
states, the attraction in the DD* system is considered to 
be about 2.7 times as strong as that of the B^B* system 
because the reduced mass of the DD* system, 967MeV, 
is about 2.7 times as light as that of the B^B* system, 
2651MeV. On the other hand, the interaction between 
the D and D* mesons is probably about the same size 
as that between the B^ and B* mesons. We argue that 
the extra attraction required for the A (3872) comes, at 
least mainly, from its coupling to the cc core, which is 
absent in the Z b systems. 

It is also an important issue that whether the charged 
partner of the A (3872) exists as a measurable peak or 
not. BABAR has searched such a state in the A — > 
it~tt J/iJj channel and found no signal [6(|. The above 
mechanism, where the coupling to the cc core plays a 
major role to bound the neutral A, is consistent with 
the absence of the charged X. 

The A (3872) exists above the open charm threshold, 
DD. Below this threshold, the cc mass spectrum is well 
predicted by a simple quark model. The model, how- 
ever, failed to predict the masses above the open charm 
threshold, especially for the DD S-wave sector. In this 
work, we also look into the cc to DD* transfer strength 
and clarify that the cc peak can actually disappear by 
the cc-DD* coupling. 

In this article, we present a hybrid picture where 
A(3872) is J PC = 1++ and consists of D°D*°, D+D*-, 
and the 2 3 Pi cc core, which stands for the x c i(2P) if 
observed. A separable DD* interaction is introduced, 
whose strength is determined so as to give a zero-energy 
bound state when it applied to the B^B* systems. The 
rest of the required attraction to form the A(3872) are 
assumed to come from the cc-DD* coupling. The cou- 
pling strength is determined so as to give the observed 
A (3872) mass. The cc core mass is taken from the quark 
model result, and the cutoff is chosen by considering the 
cc core size. As we will discuss later, the behaviors of the 



A (3872) do not depend strongly on the detail of the inter- 
actions. Main parameters of the present model are the 
overall strength of the two-meson interaction and that 
of the coupling, which are essentially determined from 
the masses of A(3872) and Zt,'s. This simple picture, 
however, is found to be consistent with many of the ex- 
periments, such as the isospin symmetry breaking in the 
A (3872) decay, the production rate of A (3872) in the pp 
collision, and the absence of Xci(2P) peak or the charged 
A, in addition to the mass of A (3872) and Z^s, which 
are the inputs. 

Let us briefly mention features of our work in compar- 
ison to those that also employ the charmonium-molecule 
hybrid model [29- 3l| . In Ref. [Hj], the hybrid structure 
of the X(3872) has been studied in the QCD sum rule 
approach by considering a mixed charmonium-molecular 
current. They found a very deeply bound A(3872), 
mx = (3.77 ± 0.18) GeV, 97% of whose component is 
a charmonium. As we shall show in SecJIIl the struc- 
ture of the A (3872) certainly depends on the binding 
energy strongly. The effective hadronic models with the 
charmonium-DZ?* molecule transition interaction have 
been used in Ref. [13, l3lll as well as the present work. 
Danilkin and Simonov [30] have studied the DD* pro- 
duction spectrum. They obtained the strength of the 
cc-DD* coupling from the heavy quarkonium decay cal- 
culated by a quark model with a small adjustment. They 
certainly found a steep rise near the threshold. We have 
studied the B -!• A(3872) K or DD* K weak decay spec- 
trum at almost same time independently in the very sim- 
ilar approach in [gj, [(32| . Here, we also introduce the in- 
teraction between D and D* , and look into the features 
of the shallowly bound X(3872). The work in Ref. [U 
they examined the effects of the Okubo-Zweig-Iizuka for- 
bidden p°J/ip and ujJ/ip channels. Since the main decay 
modes of the A(3872) are the A(3872) n+ir-J/tp and 
A(3872) — > 7T + 7r~7r° J/tjj, their inclusion is certainly im- 
portant. As shown in [31] . however, this effect on the 
pole position seems rather small. To avoid the compli- 
cation wc discuss it elsewhere. In the present study, we 
have introduced the attractive interaction between D and 
D* mesons with the coupling strength being consistent 
with the observed Z{,(10610) and Z&(10650) masses. This 
point is new to the previous two studies and we con- 
sider that we can successfully draw the consistent picture 
of the observed exotic hadrons A(3872), Zf,(10610) and 
Z fc (10650). 

One of the authors (S. T.) has studied the A(3872) 
using a quark potential model by introducing an extra 
(qq) pair to a cc system [63j and found a shallow bound 
state of qqcc with J PC = 1 ++ . Recently, an elaborate 
study has been done in the quark potential model [64| . 
They have performed the coupled channel calculations in- 
cluding two and four-quark configurations using the 3 P 
model and found a good agreement with the experimen- 
tal data. The purpose of the present work is to make 
the situation of the A (3872) clearer by studying the role 
of the cc core state, which couples to the D°D*° and 
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D + D*~ molecular states, with a simple hadronic model. 
This approach will complement the picture given by the 
quark model approach. 

This paper is organized as follows. In Scc|TTl the cal- 
culation of the X(3872) state is given. In Sec lIIIl we 
discuss the transition strength of the weak decay of B 
meson: B -> X(3872)A' or DD* K using the Green's 
function approach. In Sec lIVl we study the effects of the 
interaction between the D and D* mesons. Finally, SecfVl 
is devoted to summary of this paper. 

II. X(3872) 

We argue that the X(3872) state is a superposition of 
the cc core state, the D a D*° hadronic molecular state, 
and the D + D*~ hadronic molecular states. So, the wave 
function of the X(3872) in the center of the mass frame 
is represented by 

\X) = a |cc) + c 2 \D°D*°) + c 3 \D+D*-) . (3) 

The D°D*° and D + D*~ molecular states are given by 

\D°D*°) = J d 3 q<p (<l)\D D*°(q)) , (4) 

\D+D*-) = J d 3 qv+(q)\D+D*-(q)) , (5) 

where q represents the relative momentum of the D and 
D* mesons. The normalization of the states are 

(D D*°(q')\D D*°(q)) = (D+ D*~ {q')\D+ D*~ (q)) 



Here (po (q) and Lp+ (q) are the momentum representation 
of the wave functions of the D°D*° and D + D*~ hadronic 
molecular states, respectively. The charge conjugation is 
assumed to be positive throughout this paper. 

We introduce a coupling between the cc core state and 
the DD* states in the isospin symmetric manner. Since 
we are looking into the low energy region, the results do 
not depend much on the shape of the interaction. Thus, 
we take a monopole-type coupling as: 



(D°D*°(q)\V\ce) = (D+ D*~ (q)\V\cc) 

JL ( A2 ^ (7) 
y/A W + A 2 ) " 



The interaction we have introduced above causes effec- 
tively an attraction for the X (3872) because its energy is 
lower than the mass of the cc core, m cS . In this section, 
we ignore the direct interactions between the D and D* 
mesons; as we will discuss later in Sec lIVl the coupling to 
the cc core seems more important to make the X(3872) 
than the direct DD* attraction. 

Here we consider only the relative S'-wave states of 
these two mesons in the non-relativistic scheme because 
the X(3872) is close to the threshold. The Schrodinger 
equation to solve is 



/ TO c g - E 



V 



V m D o + %,o + -2— 

V 



V 

E 

m D + + m D ,- + 



2/H 



E, 



Ci |cc) 
c 2 \D°D*°) 
K c 3 \D+D*- 



'0 N 
= 10 



(8) 



with 
1 



1 



1 



1 



1 



1 



(9) 

/io m D o ni£)*o fi + rn D + m D ,- 

Since the interaction we employ is separable, we can 
solve this Schrodinger equation analytically. The bound 
state energy is obtained by solving the following equa- 
tion. 



with 
Fo(E) 



m cs - E - F (E) - F+{E) = 0. 



d 3 q 



(10) 



,A 3 / 2 



( mD0+mD , Q + jLy E w+^) 

(11) 



and 
F+{E) 



d 3 q 



,A 3 / 2 



For later convenience, we define ao and a+ by 



2 

2 no 



and 



= m D + + m D — — mx 



(12) 



(13) 



(14) 



where mx represents the observed mass of the X(3872). 
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In order to obtain the numerical results, we use the 
D meson masses given in the 2012 Review of Particle 
Physics [|| (TableO}. Since we have introduced the isospin 
symmetric interaction V in Eq. ([7]), the only origin of 
the isospin violation in the present model is the mass 
difference between the charged and neutral D and D* 
mesons. 

As for the cc core state, we consider that it corre- 
sponds to the J PC = 1 ++ charmonium state with the 
mass m c5 = 3.950 GeV, the closest cc core to A. This 
value is taken from the Godfrey and Isgur's results of the 
quark potential model calculation for the 2 3 Pi cc state 

In the following, we will show the results as well as 
their dependence on the various assumptions. 

There are two free parameters in the present model: 
the cutoff A and the dimensionless coupling constant g. 
We take typical hadron sizes for A: e.g., A = 0.3 GeV, 
0.5 GeV and 1.0 GeV. Then, for a given A, the coupling 
constant g is determined so that the model reproduces 
the observed mass of the A(3872), namely, 3.87168 GeV. 
The results are given in Table |TTJ The wave function we 
have obtained: 

\X) =a \cc) + c 2 l^* ) + c 3 \D+D*-) 

=ci |cc) + c /=0 \DD*;I = 0) + c J= i \DD*;I = 1) . 

(15) 

The values of c's are shown in Table IIIII for each of the 
cutoff values. It seems that the overall feature of the ad- 
mixture of each component does not depend much on the 
value of A, which is not surprising because a very shal- 
low state does not depend much on detail of the potential. 
The main component of the A(3872) state is \D°D*°), 
reflecting the fact that the mass of the A(3872) is only 
0.16 MeV below the D°D*° threshold. The amplitude 
of the \D + D*~) component is much smaller because the 
D+D*- threshold is 8.22 MeV above the mass of the 
A(3872). The size of the isospin symmetry breaking we 
have obtained for the averaged mass, 3.87168 GeV, seems 
to be roughly consistent with the one estimated from the 
experiments given by Eq. ([2]). 

Let us emphasize that we have obtained a sizable 
amount of the |cc) component. In the present scheme, 
the origin of the attraction is the coupling between the 
DD* component and cc core state. So, it is natural to 
have a rather large |cc) component in the A(3872) state. 

The mass of the cc core is taken from the quark model 
calculation. It may differ because of the model assump- 
tion. The slight change of its mass, however, does not 



change our results much. For example, for m cE — 3.93 
or 3.97 GeV, which we take because ± 20 MeV is typ- 
ical ambiguities of the quark models, the cc component 
becomes 11% and 7%, respectively. There is no drastic 
change in the results. 

We take only this state as the cc core state in this pa- 
per because it has the closest mass to that of X(3872). 
The mass of the 1 3 Pi cc state is, for example, around 
3.5 GeV and therefore its coupling to the A (3872) will 
be suppressed 64| . In our calculation which includes the 
lower cc state with the same coupling size, the probabil- 
ity of the 1 3 Pi cc state is found to be about 1/20 of that 
of the 2 3 Pi cc state. The existence of another core com- 
ponent may change the nature of the 7-decay of A(3872), 
where a large cancellation occurs and results are very sen- 
sitive to the wave function [6(1 [62}. We, however, look 
into such observables elsewhere, and concentrate on the 
bulk feature of A (3872) in this work. 

Experimental uncertainty of the X(3872) mass still ex- 
ists. So, we solve the system also for rax — 3.8687 
GeV. This mass is the one determined from the neutral 
B meson decay data, and the lightest mass among the 
ones given by the experiments. Now the binding energy 
becomes 3.14 MeV instead of the one corresponding to 
the average mass, 0.16 MeV. The value of g to form the 
lighter mass becomes 0.05625, which is 1.1 times as large 
as that of the average mass, 0.05110. In order to form 
the more deeply bound A, the dimensionless coupling 
constant g is required to be larger. The coefficients of 
the wave function are listed in Table IIIII The size of the 
cc-core component also becomes larger: it changes from 
(0.293) 2 ~ 0.086 to (0.522) 2 ~ 0.272 as m x changes from 
3.87168 to 3.8687 GeV for the case of A = 0.5 GeV. The 
size of the cc-core component in the A (3872) is found 
to be sensitive to the binding energy of the state. The 
amount of the isospin symmetry breaking depends also 
on the binding energy of A; The symmetry breaking oc- 
curs because of the difference of the binding energies of 
the A (3872) from the two thresholds, i.e., D°D*° and 
D+D*-. For ra x = 3.87168 and 3.8687 GeV,_the ratios 
of the size of the isovector to the isoscalar DD* compo- 
nents are 0.315 and 0.057, respectively. When the mass 
of the A(3872) becomes smaller, namely, the binding en- 
ergy becomes larger, the effects of the threshold differ- 
ence becomes smaller, and the isospin violation becomes 
smaller. 

Let us show the shape of the obtained wave functions. 
The explicit expressions of the wave functions in the co- 



TABLE I: Meson masses and the thresholds, 
are in GeV. 



All the entries 



m D o 


m D+ 


m D ,o 


m D «- m D o+m D ,o 


m D + + m D »- 


1.86486 


1.86962 


2.00698 


2.01028 3.87184 


3.87990 



TABLE II: The values of the dimensionless coupling constant 
g for each value of the cutoff A in units of GeV. The mass of 
X(3872) is rnx = 3.87168 GeV. 



A [GeV] 


0.3 


0.5 


1.0 


9 


0.05435 


0.05110 


0.04835 



5 




0.12 




FIG. 1: The DD* components of the X(3872) wave function 
for the m x = 3.87168 GeV and A = 0.5 GeV case. D°D*° 
wave function, r<p(r)o, is plotted by the solid line, and D + D*~ 
wave function, rip(r)+, by the dashed line. 



FIG. 2: The DD* components of the X(3872) wave function 
for the m x = 8687 GeV and A = 0.5 GeV case. The legend 
is as for Fig. [T] 



ordinate space are 

'7T\l/2 AT, 



A 2 - od v ' 



(16) 



and 

np{r) + = 

with 



7T\l/2 JVj. 



(e- Q+r - e- Ar ) , (17) 



iV = 2/x, 



g 

VA 



(I) 



N+ = 2/X+ 



(18) 



Each of the neutral and the charged DD* components 
of the wave function of the bound state with A = 0.5 
GeV and m x = 3.87168 GeV is shown in Fig. [IJ It is 
also found that the radius of the D + D*~ component is 
much smaller than that of D°D*°. In Fig. [21 we show the 
wave function of the AT(3872) also for mx = 3.8687 GeV. 
One finds that the size of the bound state, especially the 
size of D°D*° component, becomes much smaller than 
that in Fig. [T]though it is still much larger than the usual 
charmonium, whose rms < 1 fm [65]. 



III. SPECTRUM 

In this section, we investigate the transition strength 
S(E) of the B meson weak decay: B — > ccK and cc 



TABLE III: Coefficients of the X wave function. 



A 


Cl 


C2 


C3 


cr=o 


C/=l 


m x 


0.3 


0.227 


-0.947 


-0.228 


-0.831 


-0.508 




0.5 


0.293 


-0.920 


-0.259 


-0.834 


-0.468 


3.87168 


1.0 


0.404 


-0.871 


-0.280 


-0.814 


-0.418 




0.5 


0.522 


-0.727 


-0.447 


-0.830 


-0.198 


3.8687 



to A" (3872) or DD*. The A" (3872) appears as a bound 
state in the spectrum. This spectrum does not corre- 
spond directly to the observed pion distribution in the 
AT(3872) — > J/ipir n experiments. By looking into the 
DD* spectrum, however, one can see that the strength 
of DD* gathers around the threshold, and that the peak 
corresponding to the cc core actually disappears. 

In this article, we assume that the observed X(3872) 
corresponds to a very shallow bound state. To have such 
a bound state, the interaction must be attractive but 
maybe a rather weak one. As mentioned in Sect. HH we 
have fixed the strength of the DD*-cc coupling, g, so 
as to reproduce the observed AT (3872) mass. In such a 
situation, the cc core state becomes a resonance appear- 
ing in the DD* continuum. Since no sharp resonance 
is observed experimentally around 3.95 GeV, the width 
of this resonance should be large. One of the issues in 
this section is whether such a 'weak' attraction can give 
a resonance with a large decay width. 

The S(E) is normalized so that the production of the 
2 3 Pi cc state by the weak decay is equal to one. The 
vertex of the weak decay process, B — ¥ cc + K, and 
the probability that the cc is in the 2 3 Pi configuration 
are factorized out. We assume that among the cc states 
produced by the weak decay, the 2 3 Pi cc state plays a 
major role to form the AT (3872) and the DD* spectrum 
up to around E ~ 4GeV because the predicted mass of 
the 2 3 Pi cc state is 3.95 GeV. Again we use the non- 
relativistic scheme with the relative 5*-wave, because the 
reduced mass of the system is about 1 GeV and we only 
consider here up to about 0.1 GeV above the threshold. 

Then, the S(E) is expressed as follows. 



S(E) = —lm(cc\G(E)\cc) 

7T 



with the Green's function; 

G(E) = 



E - H + ie 



(19) 



(20) 



Here, E represents the energy transfer and H is the full 
Hamiltonian of the cc-core and DD* system. The state 
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4.00 




3.88 



3.92 3.96 
E [GeV] 



4.00 



FIG. 3: The transition strength S(E) with the cutoff A = 0.3 
GeV and the mass of the X(3872) m x = 3.87168 GeV. The 
cc -> V(3872) strength is 0.051, The S(E) is plotted only for 
the energy above the D°D*° threshold. The contribution to 
the bound state is not plotted here because it does not have 
a width. 



|cc) represents the center of mass system of the cc state 
with the normalization (cc\cc) = 1, This normalization 
leads the energy sum rule 



J dE S(E) = 1 



(21) 



Using the free Green's functions and the interaction given 
in Eq. (O, the Green's function is represented as follows. 



G(E) = G? + GlVGlVG\ + G^VG^VGl 



G\{E) = 



E 



■ is 



(22) 



(23) 



FIG. 4: The transition strength S(E) with A = 0.5 GeV and 
mx = 3.87168 GeV. The cc ->■ X(3872) strength is 0.087. 
The legend is as for Fig. [3] 




3.92 3.96 
E [GeV} 



4.00 



FIG. 5: The transition strength S(E) with A = 0.5 GeV and 
mx = 3.8687 GeV. The cc -> V(3872) strength is 0.269. The 
legend is as for Fig. [3] 



G° 2 (E) = 



E — ni D o — m D ,a — j— + ie 



G° 3 (E) 



E 



m D + — m D — 



2M+ 



IE 



(24) 



(25) 



The calculated transition strength for the cutoff A = 
0.3 GeV with the mass of the AT(3872) m x = 3.87168 
GeV is shown in Fig. [3J The spectrum has a sharp cusp 
above the D°D*° threshold. The resonance which corre- 
sponds to the x c i(2P) becomes very broad. The bound 
AT (3872) is not plotted in the figure because it does not 
have a width in the scheme. If we consider the experimen- 
tal inaccuracy of the energy and the A' (3872) — > J/ipTrn 
decay width, the bound A' (3872) peak and the threshold 
cusp will be merged into one single peak, which corre- 
sponds to the observed AT (3872) in the J/ipmr spectrum. 
By integrating S(E) to the DD* continuum state, one 
can obtain the transfer strength from the cc-core to the 
bound state. In this case, the former is 0.949 while the 
latter is 0.051. The cc core state of the bare mass of 3.950 
GeV becomes a resonance state of E = (3.974 — |0.067) 
GeV; its peak position is by 24 MeV shifted upward. 



In Fig.|4l we show the transition strength for the cutoff 
A = 0.5 GeV. The spectrum is almost flat at around 
E = 3.95 GeV. In the case of this harder cutoff, the cc 
core state couples to the DD* continuum of more wider 
energy range. As a result, the bump around 3.95 GeV 
found for the A = 0.3 GeV case disappears. The pole 
moves to E = (3.971 - §0.147) GeV. The strength from 
the cc-core to the bound state becomes slightly larger, 
i.e., 0.087. 

Let us now show the effect of the difference in the bind- 
ing energy. In Fig.[5j we plot the transition strength S(E) 
in the case of the cutoff A = 0.5 GeV and the mass of 
the A" (3872) is m x = 3.8687 GeV, i.e., the more deeply 
binding case. The S'-wave threshold cusp becomes much 
smaller as the bound state position moves away from the 
threshold. The transfer strength to the bound state in 
this case is 0.269, much larger than the previous cases. 

Since there is no peak is found around E — 3.95 GeV 
experimentally, the A = 0.5 GeV or more is favorable 
in that sense. This corresponds to the hadron size ~0.4 
fm, which is a reasonable value. By setting cutoff of this 
size, the shallow bound state and the large decay width 
for the cc peak can be realized simultaneously. In the 
following calculation, we use A = 0.5 GeV. 
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IV. EFFECT OF THE INTERACTION 
BETWEEN D AND D* 

In this section, we introduce the interaction between 
the D and D* mesons. First let us make a rough estimate 
of the size of the DD* attraction using the information 
from the B^B* system. Each of the BB* and B*B* 
systems has a resonances by about 2~3 MeV above the 
thresholds, Zf,(10610) and Zf,(10650). Their masses are 
10.6072 and 10.6522 GeV, respectively, and their mass 
difference is 45.0 MeV. The corresponding thresholds, 
BB* and B*B*, are 10.6048 and 10.6504 GeV, respec- 
tively, and the energy difference of these two threshold 
is 45.6 MeV. The experiment strongly suggests that the 
two-meson attraction is barely strong enough to make a 
zero-energy bound state (or somewhat weaker), and that 
there are almost no mixing between the BB* and B*B* 
states. The physical origin of the two-meson interaction 
is probably the light- meson exchange and/or the gluonic 
interaction. In ether case, the strength of the two-meson 
interaction for the DD* system is about the same size, or 
at least a similar size, for the BB* system, because the 
bosons exchanging between the light quarks is considered 
to give the largest contribution. So, as the two-meson 
interaction between the D and D* , we employ the one 
with the strength which gives a zero-energy bound state 
for the BB* systems. 

We use the Yamaguchi separable potential [68| for the 
interaction, namely, 



(MM'(q)\U\MM'(p)) = ^ v 2 



A 2 



A 2 



A 2 



A 2 



(26) 



where A is the cutoff, and A is the strength of the in- 
teraction. For simplicity, we take the same value for the 
cutoff A in Eq. fl26") as that of Eq. Q in the following 
calculation. In order to give a zero-energy bound state, 
the strength should be 



A 



A 



Ma/a/' 



(27) 



with the reduced mass of the system, fi mm' • For 
A=0.5 GeV and Hbb* =2.651 GeV, this strength becomes 
0.1886, which we denote Xb below. As for the DD*, the 
required strength to have a zero-energy bound state be- 
comes 0.5712. 

First we examine the 7r-exchange interaction if it can 
be the origin of the attraction which makes the Y(3872) 
or Zb- The 7r-meson exchange between the light quarks 
can be expressed as 



V r 



OPEP 



MAP 



(r • r)(<7 • a)V 



n 2 e -^ r 



rat m^r 



(28) 



while for the interaction between the light quark and the 
anti-quark, the factor (r • r)(<7 ■ a) becomes — (r • f)(cr • 



a), because the isospin generator for the anti-quark, f, 
corresponds to the transpose of the isospin generator of 
the quark, *t, as — f = *t. The range of the interaction, 
/i, can be obtained by using the vector meson mass, my, 
and the pseudoscalar meson mass, mp, as [x 2 = m 2 — 
{my — mp) 2 . It becomes a complex value for the DD* 
system, // 2 =-(44.5MeV) 2 , whereas /i 2 =(127 MeV) 2 for 
the BB* system. The strength V Q is about 1.3 MeV \&^. 

We assume that AT(3872) is J PC =l ++ state. Suppose 
the state consists of D°D*, then the factor (— (r-f )(c-<t)) 
is +1, whereas it becomes +3 if the state is the isospin 
state. As for the Z^'s, we assume they are /(J p ) = l(l + ) 
states. Then the factor becomes also +1 for Z& (Table 
IIV[) . Both of the values corresponds to those obtained 
from the heavy meson effective lagrangianjg^, [70[. As 
was pointed out in Ref. [69| . the Yukawa term and the 
^-function term in eq. (|2"5)) tend to cancel each other. 
Empirically speaking also, the overall effects of the OPE 
interaction is expected to be small because the factor a ■ a 
will cause a mixing between the BB* and B*B* states. 
It has been reported that the OPE interaction (with the 
tensor term and higher order partial wave states) makes 
a bound state[7(|. There, however, they found that one 
bound state below the BB* threshold and one resonance 
above the B*B* threshold rather than two similar res- 
onances. The spin dependence of OPEP seems incon- 
sistent with the B^*'B* experiments, where the energy 
difference of the two peaks is almost the same as that of 
the two thresholds. 

Thus, it seems in the B^B* systems, the two- 
meson interaction is considered to be, at least mostly, 
spin-independent isospin-independent and is just strong 
enough to make the almost-zero-energy states. In sum- 
mary, the DD* interaction we employ is: 



(D°D* (q)\U\D D*°(p)) = (D+ D*~ {q)\U\D+ D*~ (p)) 
—A / A 2 \ ( A 2 
~ A 2 



A 2 



A 2 



(29) 



with A = Xb- Though, we look into the effects of the 
DD* attraction by changing the value of A from As. 

Let us go back to the Schrodinger equation, which now 
includes the two-meson interaction, U : 



TABLE IV: The spin- isospin matrix elements of the OPEP by 
the two-meson states: DD* J PC = 1 ++ and B (,) B* I(J P ) = 
1(1+). 



(-(t ■?)(*■*)) 


D u D *u D + D *- 


D+D*- 


1 2 

2 1 


(-{r ■ r){a ■ o)) 


B+B M B*+B M 


B + B*» 
B*+B*° 


1 2 

2 1 
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/ m cE - E 

V m D o 

\ V 



m D *o 



V V 

U -E 

ni D+ +m D ,- +JL + U-E, 



p 

2^0 



Ci \cc) 
c 2 \D°D*°) 
,c 3 \D+D*-\ 



(30) 



Now the model has one more parameter, A, which stands 
for the coupling strength of the interaction between the 
D and D* , in addition to the cutoff A and the coupling 
strength between the cc core and the two-meson state, 
g. The A = limit corresponds to the results in Sec. [TT1 
where we determined the coupling strength g so as to re- 
produce the observed A (3872) mass without introducing 
the direct DD* attraction. We call that value go in the 
following and use it as a reference. 

In order to study the effects of the interaction between 
the D and D*, we change the value of the coupling con- 
stant A. For a positive A, g should be smaller than go in 
order to reproduce the observed mass of the A (3872). Or, 
equivalently, when (g/go) 2 < 1, one has to take A > to 
compensate the weakened coupling. At the 5 = limit, 
the A~(3872) becomes a pure D° D*° hadronic molecu- 
lar state. There is no charmonium component nor the 
D+D*- component in the A (3872). There will be a 
similar bound state in the D + D*~ system also, and the 
cc core becomes a sharp resonance at around 3.95 GeV. 
We consider the actual situation is in-between of the two 
A = and g = limits. 

In Fig.[6j we show the size of each of the cc, the D°D*° 
and the D + D*~ components in the A (3872) wave func- 
tion in our calculation. For each values of (g/go) 2 > we 
re- adjust the value of A to fit the mass of the A (3872) to 
be 3.87168 GeV. In Fig. [71 we also plot the sizes of each of 
the isovector and the isoscalar DD* components. As the 
interaction between the D and D* becomes larger (i.e., 
(g/go) 2 becomes smaller), the isovector DD* component 

1.0-. 
0.8- 
0.6- 



0.2- 

0.0 -| — , , , 

0.0 0.2 0.4 0.6 0.8 1.0 

(g/g ) 2 

FIG. 6: Probability of each components in X(3872). We take 
the mass of the V(3872) m x = 3.87168 GeV with the cutoff 
A = 0.5 GeV. The solid line shows the size of the cc compo- 
nent in V(3872), the dotted line shows that of the D + D*~ 
and the dashed line shows that of the D°D*°. 



D D s,ar bar (l=0) 

- . . - D D SHrba, (l=l) 











' I I — I 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

(g/g ) 2 



FIG. 7: Probability of each components in V(3872). The pa- 
rameters are the same as those in Fig. [(J] The solid line shows 
the size of the cc component in V(3872), the dash double dot- 
ted line shows that of the isovector DD* and the dash dotted 
line shows that of the isoscalar DD* . 



in the A (3872) wave function becomes larger while the 
isoscalar DD* component reduces to 0.5. 

As was mentioned before, experimentally the isovector 
component seems to be about one forth of the isoscalar 
component (see eq. ©). Also, the production process 
of A(3872) suggests that there should be a sizable cc 
component. From Fig.[7j one can find that these require- 
ments are fulfilled when (g/go) 2 is close to 1, namely the 
A = limit. 

When the DD* interaction is switched on, and its 
strength becomes A = As, the coupling to the cc core 



0.8- 




(g/g ) 2 



FIG. 8: Probability of each components in V(3872). We take 
the mass of the V(3872) m x = 3.8687 GeV with the cutoff 
A = 0.5 GeV. The legend is as for Fig. [3 
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FIG. 9: The calculated transition strength S(E) with A = 0.5 
GeV and m x = 3.87168 GeV with (g/go) 2 = 0.699. The 
cc -> X(3872) strength is 0.056. The legend is as for Fig. Q] 

becomes g = 0.0427315, which corresponds to (g/go) 2 = 
0.699. This point also gives the appropriate size of the 
isospin symmetry breaking as well as the sizable cc com- 
ponent. There each of the components of the AT(3872) 
wave function is: 

\X) = 0.237 |cc) - 0.944 \D°D* Q ) - 0228 \D + D*~) 

= 0.237 |cc) - 0.829 \DD*;I = 0) - 0.506 \DD*; 1=1). 

(31) 

This result means that about 6% of the AT (3872) is the 
charmonium, about 69% is the isoscalar DD* molecule 
and 26% is the isovector DD* molecule. Provided that 
the rhs of Eq. ([2]) corresponds faithfully to the ratio of the 
isovector to the isoscalar DD* molecular components in 
the X(3872) wave function as it is, the state expressed by 
Eq. (|31l) is consistent with the experiment. This situation 
seems to depend on the (g/go) 2 value only mildly. 

We have also solved the system where the mass mx = 
3.8687 GeV, namely, by about 3 MeV more bound case. 
The components in such a case are shown in Fig. [8] 
Here, the cc component is much larger than that of 
mx = 3.87168 GeV. The size of the cc core compo- 
nent is sensitive to the value of the binding energy. To 
make the mass mx = 3.8687 GeV, the strength becomes 
g = 0.04873, which corresponds to (g/go) 2 = 0.750. The 
core component becomes large in this situation, though 
the isovector component becomes somewhat smaller. 

In Fig. [3J we plot the transition strength S(E) for the 
A = 0.5 GeV and m x = 3.87168 GeV with (g/g ) 2 = 
0.699 case. Also when the DD* interaction is introduced, 
it is found that the strength gathers close to the thresh- 
olds. The strength to the A%3872) is 0.056. The peak 
around the cc core disappears due to the coupling be- 
tween the two-meson states and the cc core. It becomes 
a resonance of E = (3.966 - f0.091) GeV. 

Thus, we conclude that in case of the AT(3872), rather 
small amount of the interaction is coming from the direct 
interaction between the D and D* mesons and that the 
rest of the attraction is coming from the coupling to the 
cc core state. Then we have the right size of the isospin 
symmetry breaking as well as a sizable cc component, 



both of which are key features to explain the experiments. 
Also, this picture is consistent with the existence of Zt, 
resonances and absence of the charged X. 



V. CONCLUSIONS 

In this work, we have studied the structure of the 
AT (3872) as well as the transfer strength from the cc 
core to the DD* scattering state. The system consists 
of D°D*°, D+D*-, and the 2 3 Pi cc core, which stands 
for the Xci(2P) if observed. We have introduced the di- 
rect interaction between the two mesons, which is just as 
attractive as the one which makes a zcro-cncrgy bound 
state if applied to the B^B* system. Namely, we assume 
that this two-meson interaction gives the Zf,(10610) and 
Zft(10650) resonances. This interaction, however, is not 
strong enough to make a bound state in the DD* systems 
alone. In this model, the coupling between the cc core 
and the DD* two-meson state is also introduced, which 
effectively produces the attraction between the D and 
D* . We assume that this coupling provides the rest of 
the attraction required to make a bound state in the DD* 
system, X(3872). Both of the interaction and the cou- 
pling are assumed to be isospin independent. The isospin 
symmetry breaking in this model solely comes from the 
mass difference between the neutral and charged D and 
the D* mesons. 

In the obtained wave function of the AT (3872), there is 
about 6% of the cc core component. This size is consis- 
tent with a rough estimate from the AT(3872) production 
rate in the pp collision. As for the DD* components of 
the A%3872) wave function, 69% is isoscalar and 26% is 
isovector; the ratio of the isovector to the isoscalar DD* 
components is also consistent with the experiments of 
the final tt 2 to it 3 decay ratio. The present work shows 
that the structure of the X(3872) is not a simple cc nor 
a simple D°D*° bound state. It is charmonium-hadronic 
molecule hybrid, which is certainly an exotic state. 

Since the cc core cannot couple to the charged DD* 
states, such as D + D*°, the present picture can explain 
why there exists no charged partners of the A%3872). 
Also, it can explain why the 2 3 P :L cc core, or Xci(2P), is 
not found experimentally though it has been predicted by 
the quark model which gives correct mass spectrum be- 
low the open charm threshold; this core couples strongly 
to the DD* two-meson state and becomes a resonance 
with a very broad width. 

In order to confirm the present picture of the X(3872), 
we consider that the inclusion of the pJ/ip and uiJ/ij) 
channels is important because the X(3872) is mainly 
observed in the AT (3872) — > pJ/ip — > irirJ/ip and 
AT(3872) — > ujj/ip — > iririrJ/ip channels. We are now 
performing such calculations and the results will be re- 
ported soon. 

Recently, Belle Collaboration reported the results of 
the radiative decays of the AT (3872) pFTj ■ They searched 
the X(3872) — > ip'-f in B decays, but no significant signal 
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has been found. On the other hand, BABAR Collab- 
oration has reported that £>(A(3872) — > i/j'j) is almost 
3 times that of #(A(3872) -s- J/ipj) [12]. To make the 
situation clear, it is useful to calculate the radiative de- 
cays of the A (3872) in the present model including the 
charmonium structure. It is left as the future study. 
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